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Integrated optic beam combiners offer many advantages over conventional 
bullc optic implementations for astronomical imaging. To date, integrated 
optic beam combiners have only been demonstrated at operating wavelengths 
below 4 yum. Operation in mid-infrared wavelength region, however, is highly 
desirable. In this paper, a theoretical design technique based on three coupled 
waveguides is developed to achieve fully achromatic, broadband, polarization- 
insensitive, lossless beam combining. This design may make it possible to 
achieve the very deep broadband nulls needed for exoplanet searching, (c) 
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1. Introduction 

The direct imaging and characterization of Earth-like planets, especially those capable of 
supporting life, is one of the outstanding goals of modern astrophysics and science in general. 
Currently, high-angular resolution astronomical imaging can be achieved interferometrically 
by combining the wavefronts from spatially separated telescopes [1] . Interferometric imaging 
and nulling require that the light from multiple apertures be combined. Although beam 
combining can be performed using bulk optics, integrated optic (10) implementations offer a 
number of important advantages. These include spatial filtering, enhanced stability, on-chip 
fringe scanning, compactness and scalability. 10 beam combiners for astronomical imaging 
were first proposed by Kern, Malbet, Schanen-Duport, and Benech in 1996 [2]. Using silicate- 
based glass 10 devices, laboratory and on-sky stellar interferograms were demonstrated at 
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astronomical H (1.5 /im - 1.8 /iin) and K (2.0 fim - 2.4 /im) bands [3-7]. In many situations, 
especially those involving nulling and exoplanet search, operation in the infrared beyond 
3 /im, where silicate-based glasses are not transparent, is required. On-sky interferometric 
measurements have been performed in the L band (3.0 fim - 4.0 /im) using a guided-wave 
device, consisting of a two-beam fluoride glass fiber coupler, but this fiber-based technology 
is not easily scalable to multiple apertures [8]. 

We have recently developed a prototype, single-mode, integrated optic, astronomical, 
beam combiner fabricated by titanium-indiffusion in an x-cut, lithium niobate (LiNbOs) 
wafer [9-11]. The device operates in the 3.2 /im - 3.8 /xm spectral region, which lies in the 
L band, and has on-chip, electro-optically (EO) controlled, fringe scanning capabilities. Our 
results confirm that 10 devices are well suited to perform the beam combining function for 
astronomical imaging. In the infrared wavelength region, a host star is normally a million 
times brighter than the planet orbiting it, which presents major difficulties when trying to 
image the planet directly. Nulling interferometry offers the possibility to overcome this prob- 
lem by attenuating the stellar light, thus enhancing the visibility of the planet. Generally, 
achromatic phase shifting and broadband achromatic beam combining functions are required 
for deep nulling, with nulling depths of 10~^ or better over the infrared spectrum. Although 
we can image the planet by nulling at a single wavelength, there are key reasons to use a 
broad wavelength band. First, there are several key biomarkers in the infrared spectrum from 
6 /im to 18 /im. Second, the total integration time needed to detect a planet increases as the 
spectral bandwidth is reduced. 

One key component for achieving a broadband deep null is the broadband achromatic beam 
combiner. Generally, there are two basic types of planar, integrated optic beam combiners: 
reversed-Y combiners and directional couplers. The operation of a reversed-Y combiner is 
achromatic by symmetry. A reversed-Y combiner, however, when used as part of a nulling 
interferometer, suffers an inevitable 3 dB loss of signal. On the other hand, directional cou- 
pler type beam combiners can capture the entire signal, and hence are theoretically lossless 
provided that both interferometric outputs are recorded and subtracted. Unfortunately, the 
operation of the directional coupler is chromatic, and this wavelength dependence prevents 
the device from achieving deep broadband nulls. Various approaches have been suggested to 
mitigate chromaticity of 2 by 2 couplers and achieve a wavelength flattened response. These 
approaches include (i) asymmetric directional couplers, where the two constituent waveguides 
differ in width but are uniform along the direction of propagation [12], (ii) Mach-Zehnder 
type interferometers that are wavelength-insensitive [13], and (iii) tapered velocity couplers, 
where the waveguide widths are varied along the propagation direction [14]. The operation 
of the devices described by (i) and (ii) above depend on the interference between modes, 
while the tapered velocity mode coupler (iii) is meant to operate adiabatically, where there 
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is little coupling between the local normal modes. The later device is also referred to as a 
mode-evolution coupler. Unfortunately these approaches do not yield an ultra-broadband 
achromatic response nor do they preserve the initial phase difference of the input beams in 
the interference term at the output. Lack of broadband chromaticity in either power splitting 
or phase can seriously degrade the deep nulls required for exoplanets search. In this paper, 
we propose and present a theoretical design of an achromatic, broadband, polarization- 
insensitive, mode-evolution, integrated optic, beam combiner suited for space-based nulling 
interferometry. The design is based on a system of three coupled waveguides along the lines 
of [15] and [16]. The proposed device can be realized using germanium-on-silicon or germa- 
nium air-bridge waveguide structures [11]. 

In Section 2, the theory of normal modes of a three coupled waveguide system is presented. 
A design of an achromatic, polarization-insensitive, mode evolution beam combiner is pro- 
posed in Section 3. In Section 4, the derivation of the coupled mode equations that describe 
the mode coupling between the local normal modes is given, and a condition for adiabatic 
operation is presented. In Section 5, a numerical design example of the proposed achromatic 
beam combiner, based on candidate technologies for waveguide fabrication in mid-infrared 
wavelength region, is presented and the device performance is numerically evaluated. The 
results of the paper are summarized in Section 6. 

2. Normal Modes of Three Coupled Waveguides 

Before we go into the details of the design of the proposed fully- achromatic, broadband, 
mode-evolution beam combiner, we first develop representations for the normal modes of a 
symmetric three coupled waveguide system shown in the Fig. 1 [17]. The waveguides shown 
in Fig. 1 are assumed to be uniform, i.e., do not change along the direction of propagation 
z. We also assume that (this assumption is valid when the waveguides are weakly coupled) 
the total scalar (TE or TM) field,£'(x, 2), of the combined three waveguide structure can 
be approximately written as the weighted sum of the normalized waveguide modes, y\ 
$2(3^5 y), ^zis^^y)-, associated with each individual waveguide expressed in the single x, y,z 
coordinate system shown in Fig. 1. Thus, 

3 

^(a;, y, z) = ^'(x, y) exp{j(]z) = ^ ai{z)^i{x, y) (1) 

1=1 

/+00 r+oo 
/ ^fix^y) dxdy = 1, and 1=1, 2, 3. (2) 
00 J —00 

where y) is a real- valued quantity which describes the spatial profile of the mode and /3 
is the associated propagation constant. The coupled- mode equations describing this system 
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are [17] 



dai 

— — = jfj^ai + jKi2a2 + jKiaOs 
az 

da2 . 

— — = jKi2ai + jP2a2 + Jfi;23a3 
az 

— = JKigfli + JK23a2 + JPsOs 
(32 



(3) 



where /3i and /^s are mode propagation constants for outer waveguides and P2 is that of the 
center waveguide, is the couphng coefficient between waveguides i and j, and Oi, 02, and 
03 are the complex-valued mode amplitudes associated with the individual waveguide modes 
in waveguides 1, 2 and 3, respectively. We further assume that the outer waveguides are 
identical and equidistant from the center waveguide, and thus K12 = ^23 and /3i = P-^. 

Considering a longitudinally varying waveguide structure, the amplitudes of the local 
normal modes are not constant as a function of position along the waveguide. These 
longitudinally- varying local normal modes will be designated by '^q{x, y; z), y; z) and 

\l'_(x,?/; z) and expressed as 



Ei{x, y; z) = ^/(x, y; z) exp{j f3i{z)z) 

= a^l\z)<!>i{x, y- z) + a^\z)<^2{x, y; z) + al\z)<l>3{x, y; z) 



(4) 



where / = -f-, — , and ^n{x, y; z), n = 1,2,3 are the local normal modes of the n individual 
waveguides. Thus we define Ai{z) as 



Ai{z) = exp{-j(3i{z)z) 



( \ 



/ = +,-, 



(5) 



and write 



^i{x,y-z)={^^i{x,y]z) ^2{x,y]z) <l>3(x, ?/; 2) ) / 



-,-,0 



(6) 



By substituting Eq. (6) into Eq. (3), we obtain the following local normal modes of the three 
coupled waveguide structure 
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The waveguide parameter X is defined as X 



A/3 



where A/3 = (32 — f3i — K13. In 



2v^Kl2 ■ 

general, the mode profile of a single local normal mode will evolve as it propagates along 
a longitudinally-varying waveguide structure. If we consider the three local normal modes 
of the combined longitudinally-varying three-waveguide structure, it is easily verified that 
Ao{z), associated with the antisymmetric mode, remains constant along the whole structure. 
The other two symmetric local normal modes are identified by their mode order, i.e., by 
the relative magnitudes of their propagation constants, /9+(2) and f3-{z), respectively. The 
lowest order mode will be identified by the subscript +, and the third mode is identified by 
the subscript — (where /3+ > Generally, power transfer occurs between the local normal 
modes of a longitudinally-varying structure. The degree of power transfer between the local 
normal modes is a function of the overlap parameter between their spatial mode profiles, 
y; z) and y; z) where z, j G { — , 0, +}. This overlap parameter is given by [16] 



< 



^i|^|*j>=y j ^*{x,y;z)—^j{x,y;z) dxdy 



(9) 



When < \I^j|^|\E'j >= 0, as occurs for a structure that does not change along the propagation 
direction z, there is no power exchange. Similarly if the structure varies slowly enough in 
z, then -^'^j{x,y;z) will be small and < > will be approximately zero. Thus 

there will be little power exchange between the local normal modes. We propose a design 
of waveguide structures that are neither invariant along z nor vary slowly along z but for 
which < >^ and i ^ j- When there is very little power transfer between the local 

normal modes, the mode propagation is said to be adiabatic. 



3. Proposed Design 

Based on the functional forms of the local normal modes of a three coupled waveguide 
system as shown in Eq. (7), we propose to use the symmetric three coupled waveguide struc- 
ture shown in Fig. 2(a) to achieve broadband achromatic beam combining. The waveguide 
parameters are given in Table 1. 

If monochromatic input beams with complex amplitudes |ai| exp(j</)i) and {a^l exp(j (f)^) , 
respectively, are launched into the outer waveguides 1 and 3, then the excited input field 
at z = can be expressed as the linear combination of the local normal modes A+(0) and 
Ao(0) of the coupled structure at z = as indicated below 



|ai|exp(j0i 



/ 1 \ 





I«i|exp(j0i)[- 



/ 1 \ 



/ 1 \ 
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|ai| exp(j0i)[— Ao(0) + y|^+(0)] 



(10) 
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I as I exp(j03 
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ksl exp(j(/)3)[ 



-1 



/ 1 \ /l\ 

+i 

V -1 / V 1 / 

lasi exp(j03)[^Ao(O) + ^^+(0)] 



( |ai|exp(j0i) ^ 


V I«3|exp(j03) / 



[--^|ai|exp(j</)i) 



y||a3|exp(j03)]Ao(O) 



-V2 



|ai| exp(j0i 



V2 



l^sl exp(j03)]A+(O) 



(12) 



At the output of the device we wish to recover the input phase difference 0i — 03 and the full 
power, i.e., |aip + |a3p, of the incoming beams. The width of the center waveguide is kept 
fixed while the widths of the outer waveguides are equal to one another at all points along z 
but are varied along z to achieve the desired dephasing A/3(z) = [32{z)—[3i{z) — ki^. Note that 
the separation between waveguide 1 and 3 is assumed to be sufficiently large so that K13 ~ 
and thus can be neglected. The coupling term, tii2{z), and the dephasing term, A/3(2:), of 
the waveguide structure are chosen as indicated in Table 1. It will be shown in Sec. 5 that 
the waveguide transition can be made nearly adiabatic (i.e., < \E'j|^|\E'j >~ 0, i 7^ j). Thus 
the power in each local normal mode will remain constant as the mode propagates along the 
length of the structure, although the shape of each local normal mode will gradually evolve 
due to the adiabatic change of the waveguide parameters. The antisymmetric mode Aq{Q) 
and the lowest order symmetric mode A+(0) at 2; = will evolve into the following forms. 



Ao{L) 



A^{L) 





1 



exp(-j / (3o{z) dz) 
Jo 



exp(-j / P+{z) dz) 
'0 



(13) 



(14) 



6 



respectively, at the output z = L. Thus by superposition the total output at z = L becomes 



(|ai| exp(j0i) - Iflsl exp(j</)3)) exp(-j / l3o{z) dz) 



+ 
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1 \ 
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-1 / 



(|ai| exp(j(/)i) + Iflsl exp(j</)3)) exp(-j j l3+{z) dz)-^ 



1 



/0\ 



1 

Therefore, the intensities at outputs of waveguides 1,2, and 3 at z = L are given 

^1 |ai|exp(j(/)i)-|a3|exp(j(/)3) |^ = ^(|ai|^ + las]^) - ^|ai| jas] cos(</)i - (ps 

: at output ports 1 and 3 
1 



-| |ai|exp(j0i) + |a3|exp(j03 



^ = gd^il^ + I'^sH + killasl cos(0i 



: at output port 2 



(15) 



(16) 



;i7) 



From these expressions we see that the total power in all three waveguides combined is 
kiP + ksP) cind thus there is no loss in power. Furthermore, the input phase difference, 
01 — 03, is preserved in the output interference term, |ai| |a3| cos(0i — 03), which is wavelength 
and polarization independent, guaranteeing fully achromatic and polarization- insensitive op- 
eration. 



4. Coupling of Local Normal Modes 

The definition of an adiabatic waveguide transition is a transition between two waveguide 
structures that takes place in such a manner that negligible power transfer occurs between 
the normal modes as they propagate from one structure to the other. Lack of mode conversion 
facilitates device design and can lead to wavelength insensitive operation as seen in Sec. 3. 
Conventionally, the suppression of mode conversion has been achieved by slowly varying the 
waveguide structure along the direction of propagation. It was recently suggested, however, 
that a slow longitudinal variation, though sufficient, is not necessary to achieve adiabatic op- 
eration. More generally, the suppression of mode conversion is rooted in the small magnitude 
of the nonadiabatic term itself [15, 16]. Thus conventional notions of adiabaticity based on 
slowly varying changes can be abandoned, and a new and different type of adiabaticity based 
on a controlled interaction can be exploited. In this section, the coupling equations governing 
the power transfer between local normal modes and the expression for the nonadiabatic term 
will be derived. 
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4- A. Quasi-vector wave equation 

In order to calculate the spatial variation of electric field Ei{x,y,z) and magnetic field 
li{x, y,z) of the optical waveguide structure shown in Fig. 2, we will use the vector wave 
equation derived from time-harmonic form of Maxwell's equations 

(V^ + A:^n2)E + 2V((Vlnn)-E) = 0. (18) 

The time dependence of the field is of the form exp{—jut). The dielectric constant e{x, y, z) 
of a waveguide is related to its refractive index n{x,y,z) by e = eon"^, where eo is the free 
space electric permittivity, and the magnetic permeability is assumed to be its free space 
value (/i = /io) everywhere, ko = u/c = 2n/XQ is the free space wavenumber, c = l/^/JIeo is 
the vacuum speed of light and Aq is the wavelength of light in free space. 

Although it is not necessary to do so, we will adopt a scalar analysis, and thus the last 
term in Eq. (18) will be neglected. Considering quasi-TE mode (i.e., Ey ^ and Ey ^ E^) 
or quasi-TM mode (i.e., E^ ^ Ey and E^ ^ Ez), the scalar (either y-component TE or x- 
component TM) electric field of the waveguide can be written as 

Ei{x, y; z) = ^i{x, y; z) exp{j/3^'h), i = x or y (19) 

where ^i{x,y] z) is the slowly varying field profile along the propagation direction z, Pq = 
noko is the nominal propagation constant and uq is the nominal effective index of the optical 
mode. If the waveguides have no abrupt changes, then the transverse field ^i{x, y; z) will vary 
very little on the scale of a wavelength, and hence the paraxial approximation, | ^ '^^^2'^'^'' | ^ 
a#^^j^£) I ^^Yid. If we substitute Eq. (19) into Eq. (18), apply the paraxial approximation, 
and replace the vector operator by the scalar Laplacian V^, we obtain 



d _ 

oz 


1 1 




2noko 


\ dx'^ 




1 1 




2nofco 


\ dx"^ 



k,{n\x^y-z)-nl) 
dy^J 2no ^ ' 



^dz^'- 2noko \ dx^ ^ dy^)^ 2n, ^^^^ 

Next if we substitute $3; = i^a; exp(— j/3Q^'*2;) or $j, = Eyexp{—jl3^^z) back into Eq. (20) 
or (21), we obtain the following equation 

-3§^Ei{x, y; z) = Bi{z) ■ Ei{x, y- z), (22) 

where i = y for the quasi-TE mode and i = x for the quasi-TM mode. The propagation 
operator Bi{z) is given by 
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It can be easily shown that the operator Bi{z) is Hermitian. For the sake of simphcity in 
exposition, we will restrict the remaining discussion to the TE-mode and will designate the 
operator By{z) simply as B{z). The analysis for the TM-mode is similar. The three local 
normal modes \l/o(a;, y; 2;), y; 2;) and of the waveguide structure shown in 

Fig. 2 are eigenfunctions of Eq. (22), i.e., 

B{z)<l!i{x,y;z) = l3i{z)^i{x,y;z), / = +,-,0. (25) 

4.B. Coupling equation of local normal modes 

In the following analysis, the derivation of the coupling equation of the local normal modes 
will be carried out. Consider a system of three coupled waveguides with corresponding local 
normal modes \l/o(a;, 2;), \I^_|_(x, y; z) and \l/_(x, z), which are assumed to be real-valued 
(without loss of generality) since the operator B{z) appearing in Eq. (25) is Hermitian. The 
local normal modes will be normalized such that 

< '^o{x,y;z)\'^o{x,y]z) > =< ^+(x, z)|^+(x, y; 2;) > 

=< ^_(x,?/;z)|^_(x,y;2;) > 
= 1, (26) 

and the inner product < f{x, y; z)\g{x, y; z) > written in bra-ket notation is defined as 

< f{x,y;z)\g{x,y;z) >= / f* {x,y; z)g{x,y] z) dxdy. (27) 



00 J —00 



It follows from Eq. (26) and the fact that \I'o(a;, z), \l/_|_(x, 2;) and \l/_(x, 2;) are real- 
valued that 

<<iliix,y;z)\-^\^iix,y;z)>=0, 1 = 0,+,-. (28) 

Since the operator B[z) is Hermitian it also follows that the local normal modes are orthog- 
onal, that is 

<^i{x,y-z)\^k{x,y]z)>=0, l^k, {/, A;} = 0, +, -. (29) 

and therefore 

d d 

< ^i{x,y]z)\ — \^k{x,y-z) > = - < '^kix,y] z)\ — \^i{x,y] z) >, (30) 

l^k, {l,k} = 0,+,- 

Radiation modes will be neglected, and thus the total electric field, E{x, y, z), in the coupled 
waveguides can be written as the following linear combination of the local normal modes 

E{x,y,z) = aoiz)-^oix,y;z) + a+{z)'^+{x,y; z) + a_{z)^^{x,y; z). (31) 
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Substituting the above equation into the paraxial wave equation (22) 



-J-^E{x, y, z) = B{z)E{x, y, z) (32) 



yields 



oz oz 



dz dz 

.da_{z) ^ . , . , .d^.iz 
-J — 7{ W_(x,?/; z) -ja-{z 



dz dz 
= ao{z)B{z)^o{x,y; z) + a+{z)B{z)^+{x,y; z) + a4z)B{z)^I/_{x,y; z) 

= aoiz)(3o{z)'^oix, y; z) + a+{z)P+{z)^ +{x, y; z) + a_{z)P_{z)'^ -{x, y; z) (33) 

where /3o(2^); fi+{z) and P-{z) are the local propagation constants associated with the local 
normal modes \l/o(x, 2;), \l/+(x, 2) and \l/_(x, 2;). Thus 

B{z)-<^o{x, y; z) = (3o{z)mo{x, y; z) (34) 
5(z)^+(x, y- z) = /3+(^)^+(x, y; z) (35) 
B{z)^_{x, y- z) = P.{z)^.{x, y- z). (36) 

Multiplying both sides of Eq. (33) by \&+(a;, 2;), integrating over x,y and using 
Eqs. (26), (28) and (29), yields 

-J— 3ao[z) < W+(a:,?/; z)| — |Wo(x,?/; z) > 

d 

-ja^{z) < ^+{x,y;z)\ — \^.{x,y;z) >= a+{z)f3+{z). (37) 

Similarly by multiplying both sides of Eq. (33) by \l/o(x,|/; z) or '^_{x,y; z) and integrating 
over X, y yields 

-J— 3a+[z) < ■^Q[x,y;z)\ — \^+[x,y]z) > 

d 

-ja^{z) < <ilo{x,y;z)\ — \'^4x,y;z) >= aoiz)l3o{z). (38) 

■ da_{z) . , . , d , , 

-3 -jaoiz) < ^_{x,y;z)\ — \^o{x,y;z) > 

d 

-ja+{z) < ^4x,y;z)\ — \^+{x,y;z) >= a.{z)/3.{z). (39) 
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Equations (37), (38) and (39) can be written in matrix form as follows 

/ ao{z) ^ 
a+{z) 



.d_ 

dz 



-j^,^{z) /3+(^) +Ji+-{z) 





/ 


ao{z) 


\ 






a^iz) 




1 


\ 


a_{z) 


/ 



where 



d 

^ik{z)=<'^i{x,y;z)\ — \^kix,y;z)>, l^k, {/, A;} = 0, +, 



(40) 



(41) 



When there is no coupling (i.e., power exchange) between the local normal modes '^q{x, y; z), 
\l/+(x, y; z) and \l/_(x, y; z), the modes are said to evolve adiabatically. From Equation (40), it 
is clear that mode evolution will occur adiabatically when ^ik{z) = 0, I k {l,k} = 0, +, — 
for all z. It follows from Eq. (41) that mode evolution will be adiabatic when the waveguide 
structure is uniform or changes slowly along the propagation direction z since for these cases 

d 



dz 



^i{x,y;z) = or ^ 0, / = 0,+,-. 



(42) 



Equation (42), however, need not be satisfied for adiabatic operation, and the more general 
condition is given by ^ik{z) =0, I k {l,k} = 0,+,—. In future discussions, ^ik{z) will 
be referred to as the nonadiabatic terms. According to the results in Sec. 2, the three local 
normal modes of the three waveguide structure can be written as 

^+{x,y;z) = a['^\z)<!>i{x,y; z) + a^^~^\z)<l>2{x,y; z) + a[^\z)<!>3{x,y; z) 
^_(x,i/;z) = a[~\z)<l>i{x,y] z) + a^f\z)<l>2{x,y] z) + a[^\z)<^3{x,y] z) 
^o{x,y;z) = af\z)^,{x,y;z)-af\z)^s{x,y;z) (43) 

where ^i{x,y; z), I = 1,2,3 are the local normal modes of each of the three individual 
waveguides. Without loss of generality, the ^i{x,y;z), I = 1,2,3 will be assumed to be 
real-valued and normalized such that 



< <^i{x, y; z)\<^i{x, y; z) >= 1, / = 1, 2, 3. 
It follows from Eq. (44) that 

d 

< ^i{x,y]z)\ — \^i{x,y]z) >= 0, / = 1,2,3. 



(44) 



(45) 



Since the two outer waveguides are identical and symmetrically placed relative to the central 
waveguide, it also follows that 

d d 

< ^i{x,y]z)\ — \^i{x,y-z) > =< $3(x, 2;)| — |$3(x, 2) > 

d d 

< ^i{x,y;z)\ — \^3{x,y;z) > =< (^3{x,y; z)\ — \^i{x,y; z) > 



dz 



< <^2{x,y]z)\<!>i{x,y]z) > =< $2(2:, z)\<!>3{x,y] z) > 



(46) 
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Combining Eqs. (43)- (46), it may be easily verified that 



d 

io+{z) =< ^o{x,y;z)\ — \'^+{x,y;z) >= 

d 

^o-{z) =< '^o{x,y;z)\ — \<i/-{x,y;z) >= 0. 



(47) 



Consequently, 



.d_ 

^ dz 



a+{z) 



/ /3o(^) 


V 





/3+(^) 







\ 




aoiz) 


\ 






a+{z) 




/ 


v 


a_(z) 


/ 



(48) 



Thus for adiabatic operation we need only consider the coupling between the symmetric local 
normal modes \l/+(x,?/; z) and \l/_(x, z). 

4-C. Waveguide structure and nonadiabatic term 

Because the derivation of the nonadiabatic term, (z), is similar for both quasi-TE and 

-TM modes, we now will only consider the quasi-TE mode and denote the corresponding 
operator B{z) as 51+2+3(2;), 51+3(2;), and -82(2;) for the coupled three- waveguide system, 
the coupled waveguide system consisting of the individual waveguides 1 and 3 alone, and 
the individual waveguide 2 alone, respectively, with the following distribution of dielectric 
constants: 



ni+2+3(x, y; z) =nli + Anl{x, y; z) + Anl{x, y; z) 

+ A.n\{x, y; z) : coupled three- waveguide 
rii+3(x, y; z) =nli + Anl{x, y; z) + An\{x, y; z) : waveguides 1 and 3 
77,2 (x, y] z) =nli + Anl{x, y; z) : waveguide 2 alone 



An1{x,y-z) 



[x,y;z) 



n\, : inside cores i = 1, 2, 3, or otherwise. 



(49) 
(50) 
(51) 
(52) 



where n^ore and n^i are the refractive indices of the waveguide cores and cladding regions, 
respectively. The waveguides and refractive index profiles are shown in Fig. 3 for a buried 
waveguide structure. The local normal modes of the coupled system and the constituent 
waveguide modes are defined as the eigenmodes of the corresponding B{z) operator. 



^1+2+3(2;) ■ ^^(a;, y] z) = l3g{z) ■ ^^^(x, y; z) 
^1+3(2;) ■ ^i(x, y; z) = I3i+3{z) ■ ^i(x, y; z) 
B2{z) ■ ^2{x, y; z) = ^2{z) ■ ^2{x, y; z). 



(53) 
(54) 
(55) 



where the index q = {+, — } is used to distinguish the lowest (+) and the other (— ) symmetric 
local normal modes of the coupled three waveguides system (/3+ > (3-). '^i{x,y;z) is the 
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lowest order local normal mode associated with the coupled waveguide system consisting of 
waveguides 1 and 3 alone. '^2{x,y', z) is the fundamental mode associated with waveguide 2 
alone. Since the B operators are Hermitian, the eigenfunctions of 51+2+3(2;) at each fixed 
z (i.e. the local normal modes) are orthogonal and the corresponding eigenvalues are real- 
valued. Without loss of generality, the local normal modes are normalized as follows: 



< ^+{x,y]z)\^i+{x,y-z) > < '^+{x,y; z)\<iJ.{x,y; z) > 

< ^.{x,y;z)\^+{x,y;z) > < .{x,y; z)]"^ .{x,y; z) > 

It is easily verified that \[^i(x, y; z) and ^1/2(2^, z) are given by 




(56) 




^'i(x, y] z) = ^["^1(3;, y, z) + $3(x, y] z)] (57) 
m2{x,y-z) = ^2{x,y]z) (58) 

where ^i{x,y;z), ^2{x,y;z) and ^3{x,y;z) are the fundamental waveguide modes of the 
individual three waveguides. For this to be true, we need to assume weak coupling between 
waveguides 1 and 3. S{z) is the overlap integral between \l/i(x, ?/; 2) and '^2{x,y; z), i.e., 
S{z) =< \l/i(x, y; z)\'^2{x, y; z) >, and thus 

< '^i{x,y;z)\'^i{x,y;z) > < ^i{x,y; z)\'^2{x,y; z) > 

< '^2{x,y;z)\'^i{x,y;z) > < ^2{x,y; z)\'^2{x,y; z) > ^ 

(59) 

When the fundamental modes of the individual waveguides are known and the coupling 
between the individual waveguides is not very strong, the local normal modes of the coupled 
system can be approximated as a linear combination of the individual waveguide modes and 
thus ^ ^ 

' '^+{x,y;z) \ ^ f "^lix, y;z) \ _ / ci+{z) ci-{z) 
'^-{x,y;z) J \'^2{x,y;z) J \ C2+{z) C2^{z) 

Substituting Eq. (60) into the definition of the local normal modes of the coupled system, i.e., 
Eq. (53), and using Eq. (59), we obtain the following eigenvalue equation for the coefficients 
Ciq{z): 

' Bniz) Buiz) \ f c,+ iz) c,4z) 

B2l{z) B22{Z) j ■ [ C2+iz) C2^{Z) 

Ci+{z) Ci^{z) 
C2+{z) C2-{z) 



(60) 




(61) 
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where the matrix representation of the B{z) operator of the coupled system and the couphng 
coefficients between constituent waveguides are, respectively, defined as 



B2i(z) B22(Z) 



where 



+ 



1 S{z) \ / 
Siz) 1 )'[ (32{z) 

i^n{z) Kuiz) 



Kij[Z) 



K2l{z) K22{z) 

ko <^,|AiV(3-j)|^,> 
2no < ^jl^j > 



AiV(l) = Anjix, y- z)+/\nl{x, y- z) 



(62) 
(63) 

(64) 



AN{2) = Anl{x,y-z) 



(65) 



If we treat the coupled waveguide system consisting of the individual waveguides 1 and 3 alone 
(i.e., \l/i) as waveguide a and the individual waveguide 2 (i.e., \l/2) as waveguide 6, then the 
coupling coefficient, Kj^, defined by Eq. (63), is equivalent to the coupling coefficient between 
waveguides a and h provided that there is no coupling between individual waveguides 1 
and 3 of the three coupled waveguide system. By taking the inner product of Eq. (60) with 
^ q{z){q = +, — ) and making use of the normalization condition Eq. (56) and Eq. (59), the 
following additional constraint on Ciq{z) is obtained 




Ci+{z) Ci-{z) 
C2+{z) C2-{z) 




(66) 



The solution of the simultaneous equations (61) and (66) can be obtained using the method 



developed in [18] and further expanded upon by Ishikawa [16 
sumed, so that \S{z)\ <^ 1. It will also be assumed that 
introduce new parameters 6{z) and (t){z) as 

k'^^{z) 



Weak coupling will be as- 



< 1 and 



/9i-/92 

P2 



tan6'(2) 
tan<^(z) 



S{z) 
Vl^S{z? 



where 



5[,{z) 



/3i + Kii - (/32 + ^22) 



2v^r^52 



12 



1^12 + 1^21 



S Kii + K22 



2(1 -S^) 1 - 



< 1. If we 

(67) 
(68) 

(69) 
(70) 
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the solution of Eqs. (61) and (66) can be expressed as 




cos^(^l±^ _sin^(£)M£) 
cosnz) \ sin^Mz^ cos^Mzm 



(71) 



where 



/3i{z)=P{z)±^6[,{zy + ^[,{zy (72) 

Piz) = 2(1 - g2) ■^^Y:^g^{i^i2 + 1^21 + S/3i + S/32). (73) 

Using Eqs. (60) and (71), the nonadiabatic term, ,^+_(z) =< >, can be expressed 

as follows: 

lde{z) cose{z) d(l){z) 



e+-(^) 



2 dz 2cos(h(z) dz 



1 d 

<^2{x,y-z)\ — \^i{x,y-z)>. (74) 



cos (hiz) dz 



The nonadiabatic term is now expressed as a derivative of local waveguide parameters, i.e., 
coupling coefficients and propagation constants. Our goal is to minimize the nonadiabatic 
term, and thus design a nearly adiabatic device. If the center waveguide is chosen to be 
uniform along its length, then 

^*2(x,|/;^) = (75) 
and thus (since < ^^21^1^1 >= — < ^i|^|^2 >) the last term in Eq. (74) vanishes. 

5. Numerical Simulation of Achromatic Beam Combiners 

In practice, a fully adiabatic transition cannot be achieved. Therefore, there will be some 
power exchange between the local normal modes \E'_|_(x,?/; z) and \l/_(x,?/; z) as they propa- 
gate down the structure. This power coupling, however, can be minimized by careful waveg- 
uide design. The design process will be illustrated in this section by considering a germa- 
nium/silicon raised strip three coupled waveguide system. Equation (74) will be used to 
estimate the magnitude of the nonadiabatic term at each z for the design example. The 
power that couples from local normal mode \E'+(x,y;z) to ^_(x,y;z), as the local normal 
mode propagates from the beginning (i.e., z = 0) to the end (i.e., z = L) oi the device, will 
be determined by solving Eq. (48) numerically. 

5. A. Candidates of waveguide fabrication in mid-infrared region 

Both chalcogenide elements [19, 20] and group IV elements, especially silicon and germa- 
nium [21], have been studied as candidate materials for mid-IR waveguides. There are several 
promising substrate materials and fabrication methods for realizing mid-infrared waveguide 
circuits. The most promising candidates include: 
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1. chalcogenide rib waveguides on chalcogenide substrates 

2. chalcogenide rib waveguides on silicon substrates 

3. silicon or germanium rib waveguide membranes with an air bridge underneath the rib 
to serve as a lower cladding layer 

4. Ge/Si heterostructure rib waveguides with the Ge rib on top of a Si substrate 

5. hollow air core ARROW waveguides with a SiGe/Si antiresonant cladding 

The widespread use of silicon-based electronics, especially Si-CMOS technology, has be- 
come a major driving force beyond the development of silicon photonics [22]. The successful 
developments of silicon photonic devices have aroused interests in deploying silicon for mid- 
infrared applications [21]. Silicon has low loss in the spectral bands 1.2 fim - 6 fim and 24 
fim - 100 /im, while multiphonon absorption prevents its use between 6 fim and 24 fim. EO 
modulation using the free carrier plasma effect becomes more efficient at wavelength beyond 
near IR, thus enhancing the functionality of silicon optical circuits in the mid-infrared re- 
gions. Traditional substrate cladding materials, such as Si02, which are suitable at near-IR 
are too lossy to be used for device operation beyond approximately 2 fim. Suspended rib 
waveguides, as shown in Fig. 4(a), consisting of silicon membrane clad below and above by 
air have been proposed as a potential solution [21]. Another promising approach is to use a 
Ge/Si heterostructure consisting of a Ge raised strip waveguide of strain- relaxed crystal Ge 
epitaxially grown upon a silicon substrate as shown in Fig. 4(b) [21]. Crystalline Ge has low 
loss from 1.9 fim out to approximately 12.5 fim. Under certain conditions, the high refractive 
index difference between Ge and Si may permit Si to be used as the lower cladding level 
even it is lossy in the 6 fim - 12 fim band, since the penetration of the optical field into the 
cladding is limited. 

We will choose the Ge/Si heterostructure raised strip waveguide as the waveguide fab- 
rication technique to illustrate our design method. The heteroepitaxial growth of multiple 
microns of high quality germanium film on silicon by multiple steps of growth and hydrogen 
annealing has been successfully demonstrated [23]. If the silicon cladding induces too much 
loss, a Ge rib waveguide with an air bridge underneath (see Fig. 4(a)) would be a logical 
choice instead. 

5.B. Numerical design and calculation results 

The schematic plot of the achromatic mode-evolution beam combiner is shown in Fig. 5 
along with cross-section of the waveguide geometry based on Ge/Si raised strip waveguide. 
The refractive indices of the germanium strip and silicon substrate in the N band are taken 
from reference [24] and are summarized in Table 2. 
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Due to relatively large refractive index difference between germanium and silicon, the op- 
tical mode of Ge/Si raised strip waveguide is well confined inside germanium strip. Thus, the 
results obtained for the buried waveguide as described in previous sections can be similarly 
applied to the waveguide geometry based on Ge/Si raised strip waveguide. Using the 3D 
structure together with a commercial, semi- vectorial 3D mode solver (RSoft BeamPROP), 
the waveguide parameters that correspond to single-mode operation in the N band at a nom- 
inal wavelength of Aq = 10 fim for both TE and TM modes are determined. The Ge strip 
height H is chosen to be 3.5 fim and the strip width W for the center waveguide is fixed to 
be 3.5 fim as well, while the widths of the outer waveguides shown in Fig. 5(a) are varied 
from 3.58 fim at z = to 3.42 fim at z = L. Since the central waveguide is uniform with 
width W = 3.5 yum, the last term of Eq. (74) vanishes, and thus the nonadiabatic term is 
reduced to 

, idejz) cosejz) d<p{z) 

2 dz 2 COS 0(2) dz ■ ^ ^ 

Under weak coupling condition, S{z) ^ 0, and thus (f){z) ^ according to Eq. (68). Therefore 
C,+^{z) is approximately given by 

Thus, nearly adiabatic operation of the device can be achieved by minimizing where 
(see Eq. (67)) 



tane(.) = ^ (7J 



and (see Eqs. (69) and (70)) 



y f^^^ _ f^ljz) + l^lljz) - (/32(z) + K22{Z)) ^^g^ 

« M£)±^ (80) 

If the outer waveguides are identical to the center waveguide and uniform along the direction 
of propagation, it can be easily verified that 

Ku{z) = K22{z) (81) 
t^i2{z) = K2liz) (82) 

according to Eqs. (63) - (65). Thus if the variation of the width of the outer waveguides 
is much smaller than their nominal widths and the gap between the outer and the center 
waveguides, then the following relations approximately hold 

Ku{z) ^ ^22(2;) (83) 
Kuiz) ~ K2l{z) (84) 
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Therefore, Eqs. (78), (79) and (80) can be approximately written as 

. . M£) (S5) 

i^'ui^) ~ i^i2{z) (86) 
t,a„9,.) . ?^ ,87) 

A three waveguide beam combiner is designed that achieves nearly adiabatic, achromatic 
operation for the TE mode at a center wavelength of Aq = 10 fim. The waveguide parameters 
(i.e., Al3(z) and ^12(2)) from z = to z = L are varied according to Table 1 such that 
is small along the length of the device to achieve desired operation. The device's operation 
at 8 fim and 12 /xm is also numerically evaluated. 

The fundamental TE and TM modes profiles for H = W = 3.5 /im, and Aq = 10 fim 
are shown in Fig. 6. The effective refractive index, Ue/f, of TE and TM modes are 3.6017 
and 3.6454, respectively. The mode profiles are used to compute the coupling coefficient, 
K12, at different center-to-center gap spacings between outer and center waveguides for both 
TE and TM polarizations. The dephasing term, A/3 = f^(ne//2 — nefj^), is also computed 
as a function of the width of the outer waveguides. Although not shown here, the coupling 
coefficients and the dephasing terms at 8 fim and 12 /im wavelengths have also been computed 
for the purpose of device verification at other wavelengths in the N band. 

The coupling coefficient ki2{z) between the outer waveguides and the center waveguide 
is restricted to positive values, while the dephasing term A/3(z) = P2 ~ Pi can be either 
positive, or negative. In order to estimate the power transferred from the ^{x,y\ z) to 
the \l/_(x, y;^) local normal mode, we use the equivalent two-waveguide system described 
by Eq. (88). 

-J- = jPaa + jKb 
az 

£=jJ^a + jA6 (88) 

where 

a{z) = a2{z) 

b{z) = V2ai{z) = V2a3{z) (89) 

01(2;), 02(2), and 03(2) are the modal amplitude of each individual waveguide, and a{z) and 
b{z) are the equivalent modal amplitudes of the straight and curved waveguides as shown in 
Fig. 5(a), respectively. The coupling coefficient K{z) between the two waveguides in the two- 
waveguide system is related to the coupling coefficient ni2{z) between the outer and center 
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waveguides of the three-waveguide system by K{z) = \/2ki2. As indicated earher the design 
approach which will be taken is to vary A/3(z) and nuiz) along the length of the device in 
a manner consistent with Table 1 and in such a way that |^ ^ 0, where tan9(z) = 
Such a design can be realized by choosing 



2ki2(z) 



K{z) = K^axsin^z) (90) 

AP{z)=-AP^aa:COS^iz) (91) 

where K^ax is the maximum coupling coefficient (which occurs at the device center z = L/2), 
APmax is the maximum dephasing term (which occurs at the beginning and the end of the 
device, i.e., at 2; = and L), and ^{z) is a monotonically (or nearly monotonically) increasing 
function of z, with ^9(0) = and "^iL) = ir, that controls the rate of change of A/3(z) and 
K,i2{z). The waveguide parameter X{z) is equal to 2k[z) ■ '^^^ nonadiabatic term according 
to Eq. (77) can be written as 

where tan ^(2;) ^ 1/X{z). Therefore, the nonadiabatic term can be expressed in terms of the 
variation of the waveguide parameter, X{z), as follows 

The coupling between the local normal modes ^^+(2) and '$__{z) for the TE mode at Aq = 
10 /im is computed with the nonadiabatic term described in Eq. (93) for the different i^{z) 
functions listed below [15]: 

7T Z 

^^{z) = — : Linear Function (94) 
L 

71 Z '2iTX Z 

'd{z) = — 0.5 sin—— : Raised Cosine Function (95) 

TT Z '271 Z 

'&{z) = 0.426 sin : Hamming Function (96) 

TT Z ^TT Z 4:71 Z 

^{z) = — - 0.5952 sin — + 0.0476 sin — : Blackman Function (97) 
ij ij ij 

The maximum coupling coefficient, Kmaxi and the maximum dephasing term, Af^^ax, are 
chosen to be 0.0068 /im~^ and 0.012 fim~^, respectively. These values are selected such that 
the single mode condition for each individual waveguide and the weak coupling condition 
between outer and center waveguides are maintained at wavelengths from 8 fim to 12 /xm, as 
the waveguide parameters from z = to z = L are varied according to Table 1. The fraction 
of power transferred (i.e., — 10 log^g P) from the \l/+(x, y;^) local normal mode to 

the \E'_(x, y;z) local normal mode is numerically computed using Eq. (48) with a+(0) = 1 
and a_(0) = 0. The result is shown for L between and 10000 /xm for the various '^{z) 
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functions and plotted in Fig. 7. Among different "dlz) functions, for L > 3000 fim, the Raised 
Cosine and the Blackman functions outperform (i.e., smaller fractional power transferred 
from y; z) to y \ z)) the others and the Blackman function is slightly better than 

the Raised Cosine function. Based on Fig. 7, we choose L = 6000 fim for our device. Better 
performance can always be achieved by using a longer value of L. For the Blackman function 
and L = 6000 /im, the coupling term, ki2{z), and the dephasing term, A/3(z), along with 
the corresponding nonadiabatic term, are shown in Fig. 8, and the power transfer 

characteristics between the local normal modes along the propagation direction z is shown in 
Fig. 9. Using computed coupling coefficients for both polarizations, together with Eqs. (90), 
(91), and (97) with Kmax = 0.0068 fim-\ A(3max = 0.012 fim-\ L = 6000 fim, and Aq = 
10 fim, the device layout is obtained. The width variation of the outer waveguides and the 
separation between the outer and center waveguides, as a function of the position z along 
the device, are plotted in Fig. 10(a) and 10(b), respectively. 

The performance of the broadband operation of the device design is also numerically 
evaluated using Eq. (48) at wavelengths of 8 fim and 12 /im and the results are summarized 
in Table 3. Based on the data, it is seen that there is little coupling between the and the 
local normal modes, and hence the device operates adiabatically. Consequently, according 
to Sec. 3, the performance of our mode-evolution beam combiner is relatively insensitive to 
the operating wavelength and the polarization, even though the coupling coefficient, ki2{z), 
and the dephasing term, AP{z), are highly wavelength-dependent. 

The sensitivity of the performance of the device due to fabrication errors has also been 
numerically investigated. Cases where small perturbations to the coupling coefficient, ^23(2), 
or the propagation constant, /33{z), introduce asymmetries, i.e., ^23(2;) 7^ ^12(2;) and ^^{z) 7^ 
(3i{z), are considered. These perturbations are modeled here by the relationships 



where the deviation 5 is a constant independent of z. After introducing these perturbations, 
the coupled mode equations (3) are solved for the output (i.e., ai(L), 02 (i^), and 03 (L)) 
assuming that the symmetric mode A+(0) = (-^ ■^)'^ (i.e., ai(0) = 02(0) = 0, 
and 03(0) = -^) was launched into the three waveguide structure. If the device operates 
adiabatically then the output should all remain in the mode, i.e., ai{L) = 0, 02 (-^) = 1? 
and asi^L) = 0. The fraction of power coupled to the other modes, i.e., 1 — 0-2 (-^) 5 is given in 
Table 4 as a function of the deviation 6. According to Table 4, the device operation is more 
sensitive to perturbations in the propagation constant than to perturbations in the coupling 
coefficient. The predictions in Table 4. however, may be a bit pessimistic in practice since 
they assume that the errors are constant along the entire length of the device. 



S 



or 6 



f3s{z) - /3,{z) 



(98) 
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6. Conclusion 

We have presented the numerical design of a broadband, polarization insensitive, achromatic 
beam combiner for the operation in the astronomical N band based on Ge/Si heterostructure 
raised strip waveguides. The beam combining is intrinsically achromatic because of the sym- 
metric arrangement of the three coupled waveguides. As opposed to a reversed-Y junction 
combiner which suffers a 3 dB loss, our device is theoretically lossless. Furthermore on-chip 
EO modulation is also possible by utilizing the free carrier plasma dispersion effect [25] in 
silicon-based waveguides. Most importantly, the technology needed to actually fabricate the 
proposed design is quite promising and plausible. We believe the realization of such beam 
combiner will increase the possibility of using integrated optic combiners for space-based 
deep nulling interferometry. 
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Table 1. Device parameters are varied adiabatically along the structure from 
z = to z = L. 



Position z 


Parameter X 


A/3 


1^12 


1^13 





— oo 
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L/2 








^max 
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+ 00 
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Table 2. Dispersion of Si and Ge at 20°C and different mid-infrared wave- 
lengths. 



Temperature 




n{Si) 


n{Ge) 




8 


3.4158 


4.0048 


20°C 


10 


3.415 


4.0025 




12 


3.4145 


4.0012 



25 



Table 3. Fractional power coupled into the TE and TM mode at L = 6000 
/xm evaluated at different wavelengths. 



Wavelength 


^_ TE Mode 


^_ TM Mode 


8 fim 


-34 dB 


-24.1 dB 


10 /im 


-82.2 dB 


-39.5 dB 


12 fim 


-33.4 dB 


-40.6 dB 
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Table 4. Fractional power coupled to other mode, i.e., 1 — a^i^L), due to im- 
perfection of the frabrication process. 



Deviation in ^23(2) 


Deviation in (^^{z) 


l-al{L) 


0.5 % 


% 


-52.0 dB 


1.0 % 


% 


-46.1 dB 


5.0 % 


% 


-32.3 dB 


% 


0.1 % 


-38.2 dB 


% 


0.5 % 


-24.3 dB 


% 


1.0 % 


-18.3 dB 
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List of Figure Captions 

Fig. 1. Three coupled channel waveguides. Waveguides 1 and 3 are identical and equidistant 
from waveguide 2. 

Fig. 2. Schematic plot of proposed achromatic beam combiner. 

Fig. 3. Cross- section view of the coupled two- waveguide system along with refractive index 
values. 

Fig. 4. Proposed candidates for mid-infrared waveguide fabrication: (a) silicon (or germa- 
nium) rib membrane waveguide (b) Ge/Si heterostructure raised strip wave-guide. 
Fig. 5. (a) Schematic of broadband achromatic beam combiner (b) Cross-section of Ge/Si 
raised strip waveguide geometry for fundamental mode calculation by beam propagation 
method. 

Fig. 6. (a) TE mode profile at Aq = 10 /im, rig// = 3.6017 (b) TM mode profile at Aq = 10 
/im, Ueff = 3.6454 with nominal design [H = W = 3.5 yum). 

Fig. 7. Power transferred from the TE local normal mode y; z) to the TE local normal 

mode I/; 2;) at Aq = 10 /im. 

Fig. 8. (Left) The variation of the coupling coefficient, ^12(2;), and the dephasing term, A[i{z), 

along propagation direction z (Right) The nonadiabatic term, (2), for the Blackman 

function for the TE mode at Aq = 10 /im. 

Fig. 9. (Top) Total fraction of power remaining when only the local normal mode is 
excited at z = 0, i.e., a+(0) = 1 and a_(0) = 0. (Bottom) Fraction of launched power in 
mode (a+(0) = 1 and a_(0) = 0) transferred to the local normal mode 
Fig. 10. (a) The width variation of the outer waveguides as a function of propagation dis- 
tance (b) The gap variation between the outer and the central waveguides as a function of 
propagation distance. 
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Fig. 1. Three coupled channel waveguides. Waveguides 1 and 3 are identical 
and equidistant from waveguide 2. 




Fig. 2. Schematic plot of proposed achromatic beam combiner. 
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Fig. 3. Cross-section view of the coupled two-waveguide system along with 
refractive index values. 



(a) 



(b) 






Fig. 4. Proposed candidates for mid-infrared waveguide fabrication: (a) silicon 
(or germanium) rib membrane waveguide (b) Ge/Si heterostructure raised 
strip wave-guide. 
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Fig. 5. (a) Schematic of broadband achromatic beam combiner(b) Cross- 
section of Ge/Si raised strip waveguide geometry for fundamental mode cal- 
culation by beam propagation method. 




Fig. 6. (a) TE mode profile at Aq = 10 /im, riefj = 3.6017 (b) TM mode profile 
at Ao = 10 /im, rig// = 3.6454 with nominal design {H = W = 3.5 /im). 
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Fig. 7. Power transferred from the TE local normal mode "^^{XjU; z) to the 
TE local normal mode \E'_(x, j/; z) at Xq = 10 /xm. 
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Fig. 8. (Left) The variation of the couphng coefficient, ^12(2;), and the de- 
phasing term, A/3 (2;), along propagation direction z (Right) The nonadiabatic 
term, (z), for the Blackman function for the TE mode at Aq = 10 /im. 
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Fig. 9. (Top) Total fraction of power remaining when only the local normal 
mode is excited at 2; = 0, i.e., a+(0) = 1 and a_(0) = 0. (Bottom) Fraction 
of launched power in mode (a+(0) = 1 and a„(0) = 0) transferred to the 
local normal mode \E'_. 
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Fig. 10. (a) The width variation of the outer waveguides as a function of 
propagation distance (b) The gap variation between the outer and the central 
waveguides as a function of propagation distance. 
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